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I. INTRODUCTION 

Integrable Hamiltonian systems play a fundamental role in the study and description of 
physical systems, due to their many interesting properties, both from the mathematical and 
physical points of view. The construction of such models represents a contribution to this 
field, and many of them have proved to be of an extraordinary physical interest. Let us 
remind here the Morse [Q] and Poschl- Teller potentials in one dimension, or the Calogero 
and Sutherland |^ potentials. These constructions have been also considered from 
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many points of view. See for instance the reviews and . 

A method to construct these systems is the use of the Marsden-Weinstein reduction 



procedure |p, or its extensions |]10[, to free Hamiltonians lying on a A^-dimensional homoge- 
neous space under a suitable Lie group. In this way (using an appropriate momentum map), 
one assures the integrability, or even the superintegrability [|rT| of the system. In the first 
case, there exist constants of motion in involution, one of them the Hamiltonian. The 
superintegrability requires more than constants of motion (not all of them in involution) 
and more than one subset of A^ constants in involution. There are good reasons to suspect 



that any integrable system may be constructed in this way, as a reduction of a free one [|I0 
so the problem to construct these systems and study their properties is a profitable and very 
interesting field. A related topic is the problem of separation of variables for the associated 
Hamilton- Jacobi (HJ) equations. As it is well known, the existence of quadratic invariants 



allows to classify and construct these systems |jT2|,|T3|, relating them in many occasions to 
subgroups of the invariance group. 

A series of articles appeared in the last years, has been devoted to the study of these su- 
perintegrable systems constructed using homogeneous spaces over the pseudounitary groups 



SU{p, q) |T^,|TB[. In particular, using the maximal Abelian subalgebras (MASA) of the corre- 
sponding algebras, one can build a family of integrable systems of arbitrary dimension, and 
present their invariants and the coordinate systems in which the HJ equation is separated. 
The explicit solutions and a unifying view of the compact Cartan subalgebra case have been 
presented in ||T6[| . 

Our aim in this article is to work in detail the low dimensional cases. The reason is 
twofold. On one side, the one- dimensional case allows an easy geometric description of the 
systems, through their phase portrait. The potentials we obtain are not new, but have been 
applied successfully in many physical models (for instance, the Poschl- Teller and Morse 
potentials). They also appear in the study of quasi-exactly solvable (QES) models p^ , p!8 



as the case of exactly solvable systems, providing examples in which, from the quantum 
point of view, the corresponding Schrodinger equation can be solved algebraically (a part 



of the spectrum for QES systems or an arbitrary number of states for the exactly solvable 
ones). On the other side, the 2-dimensional case can be studied from the point of view of 
variable separation, and we can solve the HJ equation in a wide class of coordinate systems, 



specially in the noncompact case p]9|j20[| . The results we present here (in the 2- dimensional 
case) should be considered in a local context. Considerations about global behavior, which 
will differ from the compact to the noncompact case, will not be addressed in this work. 

The article is organized as follows. In Section 2 we present a concise description of 
the method used to construct these Hamiltonian systems. Section 3 is devoted to the one- 
dimensional case, while the 2-dimensional case is studied in Section 4. In each case we present 
the list of all the conserved quantities for these systems in terms of the generators of the 
corresponding algebras, and the explicit form in the chosen coordinate system. Conclusions 
and further outlook of this research are discussed in Section 5. 

II. INTEGRABLE HAMILTONIAN SYSTEMS AND PSEUDOUNITARY 

GROUPS 

The results presented in this section are a summary of the contents of ||l4| , p!5| . We will 



include some of them to set the notations which will be used in the following sections. 
We will consider the free Hamiltonian [n, u = 0, . . . , N = p + q — 1): 

H = Acg^'p^p. (2.1) 

(the bar denoting complex conjugate) defined in the Hermitian hyperbolic space (with co- 
ordinates y'^ G C, satisfying Qfiuy^y" = 1, and conjugate momenta p^): 

SUip,q)/SUip~l,q)xUil) (2.2) 



whose geometry is described in ||2T|. The real constant c is related to the sectional curvature 
of the Hermitian space. See also [|I6| for a detailed analysis of this space and its properties. 
Using a maximal abelian subalgebra of su{p,q) |2^, we carry a reduction procedure 
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0], in order to obtain a reduced Hamiltonian (which is not free) in the reduced space, a 
homogeneous SO{p,q) space: 

H = c (Ig^-'Ps.Ps-^ + V{s)) (2.3) 



.2 

where V{s) is a potential depending on the real coordinates s'^. The set of complex coor- 
dinates is transformed in the reduction procedure into a set of ignorable variables 
(which are the parameters of the transformation associated to the MAS A of u{p, q) used in 
the reduction) and the coordinates with the constraint g^vS^s^ = 1. 

If y^, fi = 0,...,N, is a basis of the considered MASA of u{p,q), formed by pure 
imaginary matrices, the relation between old {y'^) and new coordinates (x^, s'^) is: 

y^ = B{xr^s\ B{x) = expix'^Y^) (2.4) 

which assures the ignorability of the x coordinates (the vector fields corresponding to the 
MASA are straightened out in these coordinates). The Jacobian matrix, J, is easily obtained. 
If: 

^'^ = ^ = (2-5) 

then: 



J ^ d{y,y) 
d{x, s) 



(2.6) 



The Hamiltonian calculated in the new coordinates is written as: 



H = c i^g^'p.Pu + V{s)^ , V{s) = pl{A^KAy^p, (2.7) 

where Px are the constant momenta associated to the ignorable coordinates x and K is the 
matrix defined by the metric g. 

Note that, to obtain these Hamiltonians, we need MASAs of su{p, q) of dimension = 
p + q — 1 (corresponding to MASAs of u{p,q) of dimension p + q). We also require that 
these MASAs have a representation in terms of imaginary matrices that allows to write the 



Hamiltonian in the form (p.7|) . Once we have chosen a particular MASA, we can obtain a 
set of invariants and also the corresponding coordinate systems in which the HJ equation 



separates. The MASAs of su{p,q) are classified in E^, and for low ranks are completely 



determined. The corresponding potentials have been obtained for SU{N) in |T^, for SU (2, 1) 
5D|, for SU{2,2) in []T3[ and for any SU{p,q), choosing as MASA one of the Cartan 



m 



subalgebras, in 0]. From now on, we will always use contravariant coordinates, but we 
will write the indices as subscripts in order to simplify the notation and avoid the use of 
unnecessary brackets. 



III. ONE-DIMENSIONAL HAMILTONIANS 

One-dimensional Hamiltonians are always integrable and the phase portrait gives a com- 
plete description of the allowed motions. We shall expose the main ideas in order to achieve 
a better understanding of the more complicate systems we will study in the next section. 
We have two cases: su{2) and 5^(1, 1). 



A. 5^(2) 

We will use as a basis for su{2) the operators Xi, X2, X3, which are given in the natural 
2x2 matrix representation by: 
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in the metric K = diag(l, 1). 

In the compact algebra su{2) there is only one class of MASAs, corresponding to the 
Cartan subalgebra (CC) |2^. A basis of a representative of this class of MASA is: 



-i 



(3.1) 



A basis of the corresponding MASA of u{2) can be chosen as: 
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(3.2) 



The relation between old and new coordinates is given by the matrix B(x) in (|2.4|) . Note 
that a change of basis in the corresponding MASA of u{2) changes only the parameters 
appearing in the potential: 



2/1 = sie 



(3.3) 



The Hamiltonian, following the general expression 

■1 



H = c 



is: 



2 2 

mi mi 



!/(,) = ^ + ^ 

^0 Si 



(3.4) 



with the constraint + si = 1. Parameterizing the circle S*^ in spherical coordinates: 
So = COS0, Si = sin0, we get the Hamiltonian (c = 1): 



H{<P) 



+ V(0), V(0) 



+ 



(3.5) 



2 ■ ■ ■ ■ cos^ (p sm' 

We have only one second order conserved quantity, the Hamiltonian, which is equal to 
the Casimir of the algebra, C, up to an additive constant: 



Qi — X2 + X|. 



(3.6) 



The square of the generator in the compact Cartan subalgebra, Ci = Xf, is constant after 
the reduction and C = Ci + Qi. 

The specific values of the real positive constants mo, mi play no essential role in the 
qualitative description of the orbits and trajectories of this system. The potential has 
singularities (in the generic case) in = 0, 7r/2, tt, 37r/2. When mo or mi are equal to zero 
we have only two singularities in 0, tt or 7r/2,37r/2, respectively. 

The particles are confined inside a sector, and there, the motion is periodic, with an 
equilibrium point (a center in the phase space) corresponding to the unique minimum of 



the potential, in tan0 = ymi/mo. The solution can be easily computed, using Hamilton 
equations. The potential is bounded from below, and the energy is always positive [E > 



(mo + rrii)'^). Though the use of HJ equation is not necessary in this context of one- 
dimensional systems, we will write down the equation in order to compare with the two 
dimensional cases. In fact, when we will make separation of variables there, we will find 
again this equation: 



1 /as' 

2 [d6 



+ 



cos^ 



+ 



E 



sm 



(3.7) 



with solution iu = cos^ > 



u 



^(b+ \jh'^ -AmlE cos 2V2Et^ 



(3.8) 



and b = rn^ — m\ + E. This solution is obviously much simpler to find if we consider the 
equation of orbits in the phase portrait of this system. 
For instance, if mo = 0, mi = 1, the solution is: 



So = cos ( 



1 



cos V2Et 



(3.9) 



that is, a system with similar solutions to a harmonic oscillator, but now the frequency 
depends on the energy. 



B. sii(l,l) 

The noncompact algebra su(l, 1) has three nonconjugate classes of MAS As, compact 
Cartan subalgebra, noncompact Cartan subalgebra and a class of nilpotent maximal abelian 
subalgebras, (MANS) |2|]. We will fix the metric to be: 

(i A 

-1 



K 



(3.10) 



and the basis {Xi,X2,X3} is given in the 2x2 matrix representation through the corre- 
spondence: 
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/. Compact Cartan subalgebra (CC) 

We choose, as representative of this class, the same matrices as in ( |3.2| ). Hence, the old 
and new coordinates are related in the same way they did in the su{2) case ( ^.31 ), and the 
Hamiltonian is now: 



H = c 



,2 ^2 



V{s) = ^-^ (3.11) 



with the constraint s\ — s\ = 1. This hyperbola can be described with a coordinate (j) 
varying in the real line: sq = cosh0, si = sinh0, and the Hamiltonian in these coordinates 
is (c = -1): 

1 2 2 

IVY) TTl 

H{<P) = 7;pI + V{<P), V{<P) = + ^3.12) 

2 cosh smh 

The second order invariant (the Hamiltonian) is: 

Qi = Xi + Xl (3.13) 

and the trivial constant associated to the MASA is Ci = Xf. Hence, the Casimir in terms 
of these two quantities is C = Ci — Qi. 
The HJ equation is: 

2\d<pl cosh2 + sinh2 ^^-'^^ 

and the solution depends on the values of E and the parameters. 

Considering different values of the parameters mo,mi we obtain three different systems. 

a) If mi 7^ the potential has a singularity in = 0. It is easy to check that, if mi > mo, 
there are no minima for the potential and all the motions are unbounded (with a turning 
point). The energy is always positive. The parameters mo, mi do not modify qualitatively 
the phase portrait or the form of the solutions. The solution can be written as {u = cosh^ 0): 

u = ^ (^-6 + ^h^ + AmlE cosh 2V2Et^ (3.15) 
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and b = ml — mf — E. If mo = 0, mi = 1, the solution is: 



So(t) = W 1 + — cosh \/2Et. 



(3.16) 



b) If mo > mi > 0, the potential has two minima, symmetric respect to the origin where 
it has the singularity. The energy is bounded from below: E > —(mo — mi)^, the value 
E = — (mQ — mi)^ corresponding to the equilibrium solution in the center of the phase space. 
The other solutions are easily calculated [b has the same value as in case a)): 
i) -(mo - mi)^ < E <0 



So 



62 + 4Emg cos2J2\E\t 



1/2 



(3.17) 



ii) E = 



So 



mn 



+ 2{ml - ml)t^ 



1 1 
mo — mf 



1/2 



(3.18) 



When > we get the solution ( p.l5[ ). 

c) If mi = there is no singularity in the potential, which has a minimum in = 0, 
with periodic motions of negative energy and unbounded motions of positive or zero energy. 
The multiplicative constant mo plays no essential role for the qualitative description of the 
system. The solutions can be read off from case b) with mi = 0. 

//. Noncompact Cartan subalgebra (NC) 

A representative subalgebra of this class has the basis (in the metric ( p.lO| )): 

-i 



\ 



(3.19) 



and we will add the matrix Yq = il to get a MASA of u{l, 1). The new and old coordinates 
are related in a slightly more complicated way: 



Uo = e*^°(so coshxi + isi sinhxi), yi = e*^°(— zsq sinh xi + Si coshxi^ 



(3.20) 



and the Hamiltonian is written in the new coordinates as: 



and the constraint s^ — s\ = 1. Using again the (f) coordinate as in the previous case, we get: 



, 1 o ,,,,, mn — m? + 2momi sinh 20 
H(u) = -Pi + Vi^), V(^) = -^^-^ 1 (3.22) 



The Casimir is written as C = Qi — Ci where Ci — X^, the square of the generator of 
the noncompact Cartan subalgebra, and 

Qi = Xl - Xl (3.23) 

which is equal to the Hamiltonian. 

We will also write down HJ equation for future references: 

1 f dS\^ ttIq — + 2momi sinh 20 /o o/in 

^^^4> ^ ^ ^ ^ 

If mo or mi are equal to 0, we obtain similar results to those of the case of compact 
Cartan MASA, an attractive potential if ml > ml and a repulsive one in the opposite case. 
If mgmi 7^ 0, the potential is qualitatively the same for all values of mo and mi. 

The potential V{(f)) has two extrema in sinh 20 = mi /mo, —mo/ mi. The first point 
corresponds to a minimum (a center in the phase portrait), and the potential takes the 
value V = — mQ. The second point is a maximum (a saddle point in the phase portrait), and 
V — ml there. The energy is bounded from below {E > —m^) and the explicit solutions are 
{u — sinh 20): 
i) -ml < E <0 

1 



^ \E\ 



momi + J{E + ml){mj - E) cos 2J2\E\t 



(3.25) 



ii) ^ = 



u 



■ ""Q '"' +4momitl (3.26) 
2momi 



iii) Q < E <m\ 
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u — — 



1 

E 



momi + J{E + ml) {mj - E) cosh 2J2\E\t 



(3.27) 



iv) E — ml 



u — —mo + e 



2V2mit 



(3.28) 



\) E > m\ 



-momi + \I{E + ml) (E - mf ) sinh 2V2Et 



(3.29) 



///. Nilpotent subalgebra (NIL) 

Though the simplest representative of this class of subalgebras is obtained in the skew- 
diagonal metric, we will use again the diagonal one, because in this way, the kinetic term is 
also diagonal. We will take as a basis: 



( . .\ 



\ 



-I —I 



(3.30) 



which is a nilpotent matrix. As in the noncompact case we will also use Yq — UIo complete 
the basis of a m(1, 1) MASA. Old and new coordinates satisfy: 



yo = e'^°((l + ixi)sQ + ixisi), yi = e'''°{-ixiSQ + (1 - ixi)si). 



(3.31) 



The Hamiltonian is: 
1 



{pI-pI)+V{s) 



V{s) 



2momi 



(3.32) 



(So + Sl)2 (So + Si)4 

with the constraint (which is the same for all the subalgebras in the su{l, 1) case, as we are 
using the same metric): — s\ = 1. The expression of the Hamiltonian in terms of the 
coordinate (c = —1) is: 

1 



H{<P) = ^pI + V^(0), V{(P) = m\e~^'^ - 2m^m^e-'"^ . 



(3.33) 



The Hamiltonian in terms of the second order operators in the enveloping algebra is 
again ({X^, Xj\ = XiXj + XjXi): 
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Qi = 2X2 + {Xi,X3}-X2 



(3.34) 



and the trivial constant Ci is equal to [Xi + X^^, with C = Qi — Ci. 
The HJ equation is: 

2 \ 90 / 



(3.35) 



We will assume uii 7^ to obtain nontrivial results. Depending on the constants ttt-q, tui 
we get essentially two classes of systems: 

1) If mQrrii < there is no extremum for the potential, and the energy is always positive. 
The solutions, with a unique turning point, are given by: 



mi 



—mo + \/E + ml cosh V2Et 



(3.36) 



2) If momi > the potential has a minimum, in = (1/2) log(mi/mo) and the energy 
is bounded from below, E > —m^. As in the first case, the values of the parameters are not 
essential if they satisfy the constraints. The solutions for the different values of the energy 
are: 



i) -ml < E <0 



\E\ 



mo + J E + ml cos v 2Et 



(3.37) 



ii) E = 



^24> ^ 4momit^ 
2mn 



(3.38) 



If > the solution is the same as (|3.36|) . 

This case completes the set of one-dimensional Hamiltonians obtained through a reduc- 
tion procedure out of free Hamiltonians, invariant under SU{p,q), p + q = 2, and defined 
over a homogeneous space of the corresponding group. In Table I, we present a summary of 
these Hamiltonians in the one dimensional case. 

In the next section we will treat the 2-dimensional Hamiltonians associated to the rank 
2 algebras su{3) and su{2, 1). 
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IV. THE 2-DIMENSIONAL CASE 



There are two pseudounitary algebras to be used to construct superintegrable Hamilto- 
nians of dimension 2, su{3) and sm(2, 1). We will treat separately both cases. 



A. su{3) 



We will use as a basis for su{3) the operators {Xi, . . . , Xs} which are given in the 3x3 
matrix representation by: 



X, 



X. 



Xv 



-i 


^0 i 
i 




^ ^ 



1 
0-10 



,X2 



X. 



^ 0^ 



i 
^0 -i 

'^0 1^ 

-10 

^ 0^ 
i 
i 



,^3 



'^0 10^ 



-10 






1 



when the metric is: K — diag(l, 1, 1). 

In the compact case there is only one MASA, the Cartan subalgebra, generated by the 
matrices: 



— z 



V 







V 



(4.1) 



and we shall use the following basis for the corresponding MASA in u{3): 
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( ■ \ 



V 0/ 



Y2 



(4.2) 



The coordinates s are related to the coordinates y in the same way as in su{2), 



(4.3) 



and the Hamihonian has also the same form of all cases using compact Cartan subalgebras 

H=\(pI+p\+ pt) + V{s\ V{s) = 4 + 4 + 4 (4-4) 

with the constraint Sq + + = 1. 

In the one-dimensional case, there was only one invariant, which was the Hamiltonian. 
In this case, we can construct three invariants, only two of them in involution at the same 



time (one of them the Hamiltonian). The system is superintegrable in the sense of |]T2 
These invariants are ||14|| : 



-Roi 


= (soPi - 


SlPof + 


R02 


= {S0P2 - 


S2Pof + 


R12 


= {S1P2 - 


S2Plf + 



■So 

S2 
I — 

■So 

■Si 



■Si/ 
S2J 
82/ 



(4.5) 



The sum of these three invariants is the Hamiltonian up to an additive constant. In order 
to study the solutions of this problem we need construct a coordinate system in which the 
corresponding HJ equation separates into a system of ordinary differential equations. As in 
||T6| we will use spherical coordinates [|l9l, defined by: 



So = cos 02 cos 01, 



Si = COS02 sm(/)i, S2 = sm(/)2- 



(4.6) 



and the Hamiltonian is written as ic 
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V(M = 4-r f ^ + ^) + ^ (4.7) 

cos^ 02 \cos^0i sm 01 y sm 02 

where the constants mo,mi,m2 are chosen to be nonnegative. 

The second order conserved quantities ( |4.5| ) (we will follow the notation Q for these 
operators) can be written in terms of the basis {Xi, . . . , Xg}: 

gi = X| + X2, Q2 = X',+Xl Q3=X^,+Xl (4.8) 

with commutation relations (the commutator is a third order element which plays no essen- 
tial role in the method): 

[Qi,Q2] = [Q2,Q3] = [Q3,Qi] 

The Casimir is: 

C = 4Ci + 2C2 + 4C3 + 3Qi + 3Q2 + SQg 

where 

Ci = Xf, C2 = {Xi,X2}, C3 = X| 

are the second order operators in the enveloping algebra of the compact Cartan subalgebra. 
The Hamiltonian is 

H = Qi + Q2 + Q3 + constant 



where Qi is the expression of Qi in spherical coordinates [IS 



^ 1 2 mn fnf 
2 cos"^ 01 sm 01 

+ ]^P<t>iP<t>2 sin 201 tan 02 

„ 2,A2 2, m?\ .2,/l2 "^2 

Q3 = tan 02 [^-p,^ + j + (2^.. + 

-^P0iP02 sin 201 tan 02 



The HJ equation is: 
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1 f dS 



+ 



rrin 



1 f dS 



+ 



m, 







+ 



E 



sm 



cos' 



cos' (pi sm (pi 



(4.9) 

and separates into two ordinary differential equations using S{(f)i, = Si{(f)i) + S2{(t>2) — Et: 

(4.10) 
(4.11) 



1 /asi 

2 U(/)i 



+ 



+ 



cos^ 01 ' sin^ 01 



ai, 



1 /a^ 

2 



+ 



+ 



sm 



cos' 



02 



where a2 = E and ai are the separation constants (which are positive). These equa- 
tions have the same form as those in (p.7|) The solutions are easily computed and can 



be found as particular cases in |T^. The potential has singularities along the coordinate 
lines: 0i = 0, 7r/2, tt, 37r/2, and 02 = 7r/2,37r/2 in the generic case. It has a unique 
minimum inside each regularity domain. An analysis of the associated dynamical sys- 
tem (Hamilton equations) shows that all the orbits in a neighborhood of the critical point 
(center) are closed and hence, the corresponding trajectories are periodic (a direct conse- 
quence of the correspondence between extrema of the potential and critical points of the 
phase space). Let us restrict to the domain < 0i,02 < vr/2, where the minimum is in 
tan 01 = yJmi/mQ, tan 02 = \Jfn2/ {rriQ + mi). The value of the potential at this point is 
(mo + mi + m2)^, hence the energy E is bounded from below {E > (mo + mi + m2)^). The 
explicit solutions are: 



cos 



cos 01 



1 

2E 
1 

2^1 



^2 + \/&2 - 4«i^ cos 2V2Et 



(4.12) 



cos' 



b1 — Aaiml 



hi - AaiE 



(^(62 cos^ 02 — 2ai) sin 2v^2a7/3i 
+2^/al[{h2 — E cos^ 02) cos^ 02 — aiY^"^ cos 2-\/2ai/3i 



(4.13) 



where 61 = ai + rn^ — m\ and b2 = E + ai 



rrin 



Let us remark that these results reflect essentially the case su{2). In fact all systems we 
can construct using Cartan subalgebras can be described in a unified way as it was shown 
in |I§|TB|. 
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B. su(2,l) 



The basis we will use is formed by the set of operators {Xi, . . . , Xg} which are given in 
the 3x3 matrix representation by: 



^ 0^ 



Xr 



-z 


^0 i 
i 
^0 

^0 ^ 
1 
1 



.X. 



i 




X, 



'^0 1^ 




^10 

^0 ^ 
i 




1 0^ 



-10 
^000 

^ 00 





According to general results s'u(2, 1) has four MASAs, two of them Cartan sub- 
algebras (compact and noncompact), one orthogonally decomposable subalgebra, with one 
nilpotent element, and one nilpotent subalgebra. We will discuss these four cases in the fol- 
lowing. Although some of these subalgebras have a simpler expression in some skewdiagonal 
metrics, we will always use the diagonal one: 



K 



V 



(4.14) 



because, in this way, the kinetic part is always diagonal. There are nine coordinate systems 
associated to 0(2, 1) free Hamiltonians: spherical, hyperbolic, elliptic (I and II), complex el- 
liptic, horospheric, elliptic parabolic, hyperbolic parabolic and semicircular parabolic p0| , p3| . 
Not all of them will separate our systems because these are not free. However, the appro- 
priate systems have been computed in [pO| and we will use their results. 
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/. Compact Cartan suhalgehra (CC) 

The compact Cartan subalgebra has a basis formed by the same two matrices we used 
in sm(3), and the same situation happens for the corresponding MASA in m(2, 1) (^4.2| ). The 
coordinates s are also related to the coordinates y as they did in the compact case sm(3) 



However, the Hamiltonian reflects the noncompact character of sm(2, 1): 

H = c(\[pI+pI- pt) + Vis)) , V{s) = 4 + 4-4 (4-15) 

where the constraint Sq + — S2 = 1 must be satisfied. 

This Hamiltonian separates in four coordinate systems, spherical, hyperbolic and elliptic 



I and II [^. We will use spherical coordinates to discuss the explicit solution. 

So = cosh 02 cos 01, Si = cosh 02 sin 01, S2 = sinh02. (4.16) 
Choosing c = — 1, we have the Hamiltonian in these coordinates: 



cosh 02 \cos2 0i sin 0i y sinh <p2 
Due to the form of the potential the constants mo,mi,m2 can be chosen nonnegative. 
The potential is regular inside the domain O<0i<7r/2,O<02<oo. It has a saddle point: 



tan 01 = y/mi/mo, tanh02 = \Jm2/ {mo + mi), if mo + mi > m2. However, due to the 
special form of the kinetic term (which is not positive definite), it is easy to check that the 
associated dynamical system has all the orbits in a neighborhood of the critical point (which 
is also a center as in the compact case) closed and again, the corresponding trajectories are 
periodic. 

The second order operators in the enveloping algebra of this MASA are: 

Ci = X^, C2 = {Xi,X2}, C3 = X2 

The quadratic constants of motion lying in the enveloping algebra of sm(2, 1) and commuting 
with the elements in the compact Cartan subalgebra are: 



Qi — + Xl, Q2 — X^ + Xq, Q3 — X^ + Xg 



with commutation relations: 



[Qi,Q2] = [Q3,Q2] = -[Q3,Qi] 



The Casimir is written in terms of these second order operators as: 



C = 4Ci + 2C2 + 4C3 + 3Qi - 3Q2 - 3Q3 



and the Hamiltonian is: 



H = -Qi + Q2 + Q3 + constant 



where Qi are the conserved quantities in spherical coordinates: 
1 2 ^0 



2 cos^ 01 sm 01 

Q2 = tanhV2 f^pji 



ml 



-p^^ sin 01 , 

COS"^ 01 



-\p<t>iP4>2 sin 201 tanh 02 



2 I / 2 2 I 
= tanh 02 I -p^^ cos 0i + ^ 



sm ! 



61 tanh ( 



\ • 2 , A 2 ^2 ^ 

-^j+sm 01 (2^.. + ^^^^ 



The HJ equations corresponding to the Hamiltonian (|4.17|) are: 



1 (dSi 



+ 



2 2 



1 /9S2 



cos^ 01 sin^ 01 

— r"2 — — — 1"^ — ~ Ol2- 
sinh 02 cosh 02 



The first one is the same as we got in 5^(3) ( [4.10|) , ai is always positive and ^2 

The solutions depend on the values of the parameters and energy 
i) E < 



U2 



2\E\ 



-b2 + \lhl + AaiEcos2\l2\E\t 



ii) E = 



19 



U2 



«! — 1712 



+ 2{ai - ml)t^. 



(4.22) 



iii) E > 



1 



U2 



2E 



&2 + V ^2 + 4ai^ cosh 2V2Et 



(4.23) 



where U2 = cosh 02, &2 = — ai + The other equation can be solved as we did in the 
previous cases. The result is: 



Ui 



1 


hi H 


1 


\hl- 




1/2 


2ai 




1 

U2 




- AaiE 





-(62'U2 + 2ai) sin2\/2a7/9i 



+2^i[{Eu2 - h2)u2 - aiY'^ cos2v^/3i 



(4.24) 



where Ui = cos^ 0i and hi = ai + ml — 
//. Noncompact Cartan subalgebra (NC) 

There is only one noncompact Cartan subalgebra. A representative can be chosen ac- 
cording to the same criteria we used in sm(1, 1) (|3.19|) , keeping one element compact and 
the other (as in sm(1, 1)) noncompact: 



V 



—I 



—I 







V 



i 
-i 



(4.25) 



and the basis for the corresponding MAS A of m(2, 1) will be: 



/ 



Yn 



\ 



Yi 



( 



Yo 



i 
-i 



(4.26) 



The coordinates are as in su(lA): 



yi = e*^^ (si cosh X2 + is2 sinh X2), 
y2 = e*^^ {—isi sinh X2 + S2 cosh X2] 
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(4.27) 



The Hamiltonian is: 



H = c(^^{pl+pl-pl)+Vis)^ 

,2 ^2\/„2 „2^ 



si + {si + sly ^ 



where we will take rriQ > and mi, 1712 can take any value, and the coordinates satisfy the 
constraint (the same for all su{2, 1) MASAs, as we have chosen the same metric in all cases): 

Sq ~^ ^1 ^ ^2 ~ ^■ 

There are two systems of coordinates in which the associated HJ equation separates, 



hyperbolic and complex elliptic |2y]. We will use hyperbolic coordinates, defined as: 



50 = cosh 02 

51 = sinh 02 sinh 0i (4.29) 

52 = sinh 02 cosh 0i 

and the new Hamiltonian is (c = —1): 

cosh 02 smh 02 \ cosh 20i / 

Note that the potential follows the same pattern as the corresponding case in su(l, 1). It 
is regular inside the domain: — oo<0i<c)O,O<02<cx3, and has also a saddle point at: 



sinh 201 = —1712/1711, tanh02 = ^J\mi/mo\, when \mi\ < |mo|. As in the previous case, the 
associated dynamical system has a center and the trajectories in a neighborhood of it are 
periodic. 

The basis for this MASA is {2X1 + X2, Xg}, and the corresponding second order elements 

are: 

Ci = (2X1 + X2Y, C2 = {2X1 + X2, Xs}, C3 = X| 

The second order conserved quantities, commuting with 2X1 + X2 and Xg, and belonging 
to the enveloping algebra of su{2, 1) are: 
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Qi — X2 — X^, 

Q2 — xl + xl — xl — Xq, 

Qs^{Xs,X,} + {X,,Xe} 



(4.31) 



with commuting relations: 

[4,4] = [4,4], [4,4] = o 



The Casimir is written as: 



and the Hamiltonian is: 



H — Qi + Q2 + constant 



Finally, the conserved quantities are expressed in hyperbolic coordinates by: 

1 m 

Qi 



ii = 



1 2 iTii — — 2mim2 sinh 20i 



cosh^ 201 



ml 



„ 2 ""0 I 

V2 = 2^^, - ^Qgj^2 - tanh^ 02 1^2 

1 o \ cosh 201 



mi — m| — 2mim2 sinh 20i 



cosh^ 201 



^ i . , ^ , / 2 2 \ cosn zcpi 

Q3 = 2 ^-1^ 20, ^p,, + ^^^P,. j - -^^^n^p^.^^^. 

, 2. i2, • uoA (m^ - /;/^)siiiri2oi + 2/;/i///2 

+mQ tanh 02 smh 20i 2 2 

tanh 02 cosh 20i 

The Hamilton- Jacobi equation separates into two ordinary differential equations: 
1 / dSi \ ^ ml — ml — 2m\m2 sinh 20i 



cosh^ 201 



ml 



CKi 



The solutions have the same 
i) E<Q 



1 

2y902/ cosh^02 sinh^02 
form we have found before. 



(4.32) 
(4.33) 



U2 



2\E\ 



62 + -4q;i£;cos 2^j2\E\t 



(4.34) 



ii) = 



22 



U2 



+ 2{ai + ml)t'^. 



(4.35) 



iii) E >0 



U2 



1 

2E 



h + \ hl - AaiE cosh2V2Et 



(4.36) 



where U2 = sinh^ <^2, &2 = + ai + ttiq. 



The solution for the other coordinate is obtained in the same way (with the change 
ui = sinh20i). The equation for this coordinate is the same as that given in formula ( |3.24| ) 
and its solutions can be found in formulas ( p.25| - p^29l) . Due to the possible different signs 
of the energy and the constant ai, one should take care of the square roots appearing in all 
the formulas. 

///. Orthogonally decomposable subalgebra (OD) 

The orthogonally decomposable subalgebra (a representative of the class) has a basis 
formed by a compact element and a nilpotent one: 

/ 



2i 



—I 



V 



V 



-i —i 



(4.37) 



and the basis for the corresponding MASA of m(2, 1) is: 



V V 



\ 



-i —i 



(4.38) 



The coordinates have also a similar form to those in su(l, 1) ( |3.31| ): 



yi = e''''{{l + ix2)si + ix2S2), 



y2 = e'''^{-iX2Si + (1 - iX2)S2). 



(4.39) 



The Hamiltonian is: 



23 



H = c(^^{pl+pl-pl)+V{s)y 

ml m^Si — S2) 2mim2 
So {si + 82)^ {si + S2) 



V{s) = ^ - 7^^^~ + -^^^ (4.40) 



with sl + sj — si = 1. 

There are four coordinate systems associated to this subalgebra, hyperbohc, horospheric, 
elhptic parabohc and hyperbohc parabohc We wiU use again the hyperbohc ones, 



defined as in (t4.29| 



The Hamiltonian is (c = —1): 

V{h,h) = ^ -i— Ule-'*' + 2m,m2e''*') (4.41) 

cosh 02 smh 02 

The potential is regular inside the domain: — oo<0i<oo,O<02<oo, and, as it 
happens in all the 5^(2, 1) cases, has a saddle point at: 0i = (1/2) log(|m2/mi|), tanh02 = 



mi/mo|, when |mo| > |mi|,mim2 < 0. The situation is the same as in all other cases in 
su(2, 1). 

The second order operators in the enveloping algebra of the MASA under consideration 
are given by: 

Ci = (2Xi + X2)^ C2 = {2Xi + X2, X2 + Xg}, Cs = (X2 + Xs)^ 
and the quadratic constants of motion: 

Qi = X| + X4 — Xg — X|, 

Q2 = (X3 + + (X4 + X,y, (4.42) 
Q3 = X2 + 2{Xi,X2 + X8} 

satisfy the commutation relations: 

[Qi,Q2] = [Q2,Q3], [Qi,Q3] = o 
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The Casimir is given in terms of these operators by 



and the Hamiltonian is: 



C = Ci + 3C2 - 3C3 + SQi - SQ. 



H = —Qi + Q3 + constant 



We can write the conserved quantities in hyperbohc coordinates: 



m. 







Q2 = e'f^ 



(Ipj^ + ml tanh^ 0, + ^-i^ (^p^ + mje 



2"^^ cosh^(^2. 

2^-401 ' 



tanh( 



Q3 = 2^01 + "^2e 



2 , ^2^-4<^, ^ 2mim2e-2<^^ 



The HJ equation is separated into the following equations 



1 /dSi 

2 [dd^ 



m 



2g-4<^i ^ 2mim2e-2'^i 



1 fdS2 

2 U02. 



cosh 02 sinh 02 



E. 



(4.43) 
(4.44) 



Equation ( |4.44| ) is integrated using U2 = sinh^ 02. The result is the same as in the 
previous case ( |4.33| ). Equation ( |4.43| ) is solved using ui = e^"^^ (the same change we use in 



the nilpotent MASA of the su{l, 1) case), and the results are essentially the same we have 
found above (see p.35| ). 



IV. Nilpotent subalgebra (NIL) 

The nilpotent subalgebra has a basis formed by two nilpotent elements (one of order 2 
and the other of order 3): 







—i —i 



^ . 



i 
-i 



(4.45) 
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and the basis for the MASA of u{2, 1) can be obtained adding to these two matrices the 
identity times the imaginary unit. 

The new coordinates are defined through: 

yQ = e^^o(so +iX2(Si +S2)), 

yi = (^ix2So + (1 - y + ^^i^ ^1 + Y + ^^i^ ' (4-46) 
The Hamiltonian is: 



H = c{^1{pI + pI-pI)+V{s)^ 



_ 2momi + ml _ Amim2So mjiAsl - 1) 

isj + S2y isl + S2)'^ isl + S2r ^ ' 

with the constraint: + — = ^■ 

We have now two separable coordinate systems: horospheric and semicircular parabolic 
||20| . We will use now the horospheric ones, defined by: 

so = 0le^^ Sl = cosh02-^0?e'^^ S2 = sinh 02 + ^^'e^^ (4.48) 

The Hamiltonian is (c = —1): 

V{(j)i, 02) = m?e"^'^2 - e~'^^^ {ml + 2momi + 4mi0i(mi0i - m2)) . (4.49) 

The potential has no singularity in the whole plain (0i,02)- It has a saddle point at: 
01 = m2l2mi, 02 = (1/2) log(mi/mo), when m^m^ > 0. The situation is the same as in all 
other cases in 5^(2, 1). 

The second order elements in the enveloping algebra of the nilpotent subalgebra are: 

Ci = (X2 + X^)\ C2 = {X2 + Xg, X4 + Xe}, C3 = (X4 + X^f 
and the constants of motion 
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Qi 



3(^3 + X,f - 2{2X, + X2, X2 + Xg}, 



Q2 



{2X1 + X2, X4 + Xe} + 6{X4, + Xg} - 3{X7, X3 + X5}, 



(4.50) 



<53 



4X1" + 3X2' - 2{Xi, X2} + 6X| + - 3X7' - {4Xi - X2, Xg} 



+ 3{X3,X5} + 3{X4,X6} 



have the following commutation relations: 



[gi,g2] = [g3,g2], [QuQz] = ^ 



The Casimir is 



c = -3Ci - 3C3 - (9i + (93 



and the Hamiltonian.: 



H — Qi — Q^-\- constant 



Finally, the second order constant of motion are given in horospheric coordinates by the 
following expressions: 



This is the most interesting case, in the sense that the others are easily reduced to the 
cases in dimension 1, while this nilpotent subalgebra does not appear in the sit(l, 1) case. 
However, the solutions are still very similar to those found before. It is worth mentioning 
here, that all the potentials we have construct (and any potential we could construct by 
using this method) are always inverse quadratic potentials in the coordinates, and the solu- 
tions have always similar forms (though they depend on the specific characteristics of these 
potentials and the constants involved). 



Q2 



Qi 



-pI^ + 4mi0i(mi0i - 7712) 
1 1 

2</'ip5i - 2^01^02 + {ml + 2momi)(/)i 

— mie~^''^^(2mi0i — 777.2) + 4777i0^ (7771^1 — 7772) 
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The HJ equation is separated into the following equations 



1 dSi 

2 U^i 



-\-m\-\- 2momi + 4mi0i(mi0i — 771.2) = cti, 



l^\\rnle-^'i'^-a,e-^'t>^^E. 



1 
2 



(4.51) 
(4.52) 



The change U2 — 6^*^^ allows to solve the second equation, and the other one is solved 
directly. The solutions are: 
i) E <Q 

1 



U2 



2\E\ 



ai + Jaf + 4mfEcos2J2|^|t 



ii) E = 



iii) E>^ 



U2 



2E 



ai + V «i + 4m? cosh 2V2Et 



The first equation gives the value of the 0i coordinate: 

nl/2 



01 



1 



2mi 



m2 H 



ai — 2momi 
a? + AmlE 



aiU2 — 2ml) sin 2-\/2/3imi 



+2mi[(Eu2 - ai)u2 - va^^l'^ cos2V2/?imi) 



(4.53) 



(4.54) 



(4.55) 



(4.56) 



In Table II, we present a summary of these Hamiltonians in the two dimensional case. 



V. CONCLUSIONS 

We have presented in this work a complete analysis of a series of one and two dimensional 
integrable Hamiltonians, which in the 2-dimensional case are superintegrable in the sense 
described in the Introduction. Though the one-dimensional case is always an integrable 
system, let us remark the importance and applications of the potentials described in Section 
2. Regarding the two dimensional ones, we have provided them with a set of conserved 
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quantities which allows to study the HJ equations in several separable coordinate systems 
and compute in some interesting cases the explicit solutions. 

The one-dimensional Hamiltonians obtained here have been extensively studied in the 
literature from many other points of view. See for instance for a recent application 
of Morse potentials. As an example of these different approaches, all of them appear in 
the classification of quasi-exactly solvable Schrodinger operators ||l7|Jl8[| , as particular types 



of these systems corresponding to the so called exactly solvable systems. Following the 



classification in the exactly solvable potentials of Cases 1 and 2 are just the ones we 
have obtained associated to 1) and its compact and noncompact Cartan subalgebras. 
The first one ( |3.12| ) is the celebrated Poschl- Teller potential. That appearing in case 3 is 



the Morse potential ( p.33|) , which we get using the nilpotent subalgebra of su{l, 1). Finally 



the potential ( p.5|) , associated to the Cartan subalgebra of su{2) is related to the modified 



harmonic oscillators appearing in |18| as cases 4 and 5. One has to take into account in this 



case, that QES potentials, as studied in [0, are defined in the line (or half-line), and we are 



working here in a sector of (see also for a study of harmonic oscillators in a sector). 
The relation is not surprising at all if one considers that QES systems in the line are related 
to the complex Lie algebra s/(2) [|r^, and the ones we get here reflect the invariance under 
su{2) and su{l, 1), the real forms of s/(2). In the QES setting, Schrodinger operators belong 
to the enveloping algebra of a Lie algebra, while in our approach, the corresponding classical 
Hamiltonians are second order Casimirs of the algebra, and hence, they are particular cases 
(exactly solvable) of the former. 

Two prolongations of this study are now in progress. One of them is the use of contrac- 
tions in Lie algebras to obtain other Hamiltonian systems associated to different algebras, 
not necessarily semisimple. In this sense, the Hamiltonians, the conserved quantities and the 
coordinate systems can be obtained by contraction p6|j27[] . The second one is to apply this 



approach to the quantum case, considering the Schrodinger equation with these potentials 
. We also plan to study the links of this theory with QES systems and the possibility of 
considering partial integrability and partial variable separation in HJ equations p9| . 
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TABLE I. One-dimensional potentials 



Algebra 



MASA 



Kinetic term 



Potential 



su{2) 



Compact Cartan 



pI 



m. 



^ ^? 



cos^ sin 



su{l,l) 



Compact Cartan 



Noncompact Cartan 



Nilpotent 



pI 



pI 



rrin 



+ 



cosh sinh 



ml — m{ + 2momi sinh 20 



cosh 20 



TABLE IL Two-dimensional potentials 



Algebra 



MASA 



Kinetic term 



Potential 



su{3) 



CC 



4>2 



+ 



<p2 \ cos^ ^1 sin^ 01 y sin^ 02 



5-^(2, 1) 



CC 



NC 



OD 



NIL 



1 ( 2 _ ''I' l 

2 \P<I>2 cosh^ 02 



cosh 



1 / mg rn^ \ _ 

1^ 02 V cos^ 01 sin^ 01 J si 



2 sinh^ 



cosh 



^02 sinh^ 02 V 



sinh^ 02 



-rrij — 2mi m2 sinh 20i ' 
cosh^ 201 , 



i ( „2 



</>l 



2 1^02 sinh2 02 



' 1 

^ - -A, (mie-4^i+2mim2e-"'»i; 



cosh 02 sinh 02 V ^ 



2 (^02 



p2 _g-2,^.2„2 



'tf e -^'^2 -e- 2<*2 m| +2momi +4mi 0i (mi 0i -m2 ) ) 
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